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ABSTRACT 

Experiments  were  conducted  at  the  NURDC  oceanographic  tower  off  San 
Diego,  California  to  describe  the  small  scale  physical  properties  in  the 
upper  ocean  and  to   determine  their  temporal  and  spatial  interrelation- 
ships.  The  measured  parameters  included  water  particle  velocities,  tem- 
peratures, salinity,  sound  speed,  sound  speed  (phase)  modulation,  sound 
amplitude  modulation  and  surface  waves.   Simultaneous  time  series  meas- 
urements are  analyzed  for  a  twenty  minute  record  taken  at  a  depth  of 
7.3m.   The  spatial  correlation  distance  for  temperature  was  found  to  be 
about  25  cm.   Sound  modulation  of  a  60  kHz  sound  beam  underwent  a  mean 
attenuation  of  about  5%.   The  correlation  distance  of  an  effective  in- 
homogeniety  due  to  the  combined  temperature,  salinity,  and  bubble  effects 
was  found  to  be  8.8  cm.   The  change  in  the  index  of  refraction  of  the 
speed  of  sound  at  60  kHz  (a  bubble  resonant  region)  was  found  to  be  9970 
induced  by  bubbles.   The  energy  density  spectra  were  dominated  by  wave 
induced  motion.   First  order  wave  theory  accurately  predicts  orbital 
water  particle  motion  and  its  irrotationality .   High  coherence  is  found 
between  wave-induced  water  particle  motion  and  sound  modulation  and  de- 
creasing coherence  between  water  particle  motion  and  salinity,  sound  vel- 
ocity, and  temperature. 
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I.   INTRODUCTION 

A .   BACKGROUND 

The  importance  of  detailed  studies  of  sound  propagation  in  the  upper 
layer  of  the  ocean  has  become  increasingly  apparent  as  the  role  of  the 
submarine  continues  to  grow  in  strategic  and  offensive  capabilities.   Also 
significant  from  a  scientific  viewpoint  is  that  the  manner  in  which  sound 
propagates  can  be  used  to  effectively  describe  properties  of  the  medium. 
Thus  in  order  to  understand  sound  propagation,  a  detailed  knowledge  of 
the  parameters  that  affect  it  and  their  interrelationships  with  each  other 
and  sound  is  required.   Wilson  [lO]  developed  an  equation  that  related 
the  speed  of  sound  in  water  to  temperature,  salinity,  and  pressure.   Lie- 
bermann  1.4]  delved  into  the  effects  of  temperature  inhomogeneities  on 
sound  propagation  and  determined  that  inhomogeneities  do  have  an  effect 
on  the  intensity  and  refraction  of  propagated  sound  waves.   Mintzer  and 
Stone  [6]  theoretically  developed  an  equation  that  related  to  the  reduction 
in  sound  intensity  of  a  propagated  wave  to  wave  number,  correlation  dis- 
tance , distance  traveled  and  changes  in  the  refractive  index  of  the  medium. 

In  November  of  1971,  Bordy  [2],  Duchock  [3],  Rautmann  [7],  Seymour 
L8] ,  and  Smith  [9]  from  the  Naval  Postgraduate  School  took  measurements 
at  the  NURDC  tower  of  various  oceanographic  parameters  in  the  near  surface 
zone  while  propagating  sound  waves  over  prescribed  paths.   The  parameters 
measured  included  sound  velocity,  sound  amplitude  and  phase  modulation, 
water  particle  velocities,  temperature,  salinity,  and  wave  height.   The 
main  purpose  of  this  paper  is  to  extend  the  results  of  the  references  just 
cited  carrying  out  a  more  detailed  study  of  these  paramenters  and  their 
interrelationships . 
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B.   OBJECTIVES 

1.  To  determine  the  interrelationships  between  the  various  para- 
meters (temperature,  salinity,  sound  velocity,  modulation  and  wave  height) 
in  the  near  surface  ocean  zone.   Only  those  processes  with  periodicity 
between  0.5  to  50  seconds  will  be  considered. 

2.  To  compare  results  with  established  equations  and  relationships. 
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II.   EXPERIMENT 

On  21  October  1971,  Bordy,  Duchock,  Rautmann,  Seymour,  and  Smith 
from  the  Naval  Postgraduate  School  conducted  experiments  at  the  Naval 
Undersea  Research  and  Development  Center  Oceanographic  Research  Tower 
at  San  Diego,  California.   The  purpose  of  the  experiment  was  to  describe 
the  small  scale  physical  properties  in  the  upper  ocean  and  to  determine 
their  temporal  and  spatial  interrelationships.   The  tower  is  approximately 
one  mile  off  Mission  Beach  in  a  mean  water  depth  of  17  meters.   The 
structure  is  made  of  steel  and  concrete  and  built  into  the  sea  floor.   An 
enclosed  upper  deck  provides  space  for  electronic  equipment  while  an  open 
lower  deck  provides  a  platform  to  place  instruments  in  the  water.   On  the 
northern,  southern,  and  western  sides  are  located  winch  controlled  in- 
strument carriages  that  can  be  adjusted  to  any  desired  depth.  (Figure  1 
and  Figure  2) . 

The  parameters  measured  during  the  experiment  included  orthogonal 
water  particle  velocities,  an  array  of  temperatures,  salinity,  sound 
speed,  sound  speed  (phase)  modulations,  sound  amplitude  modulation  and 
surface  waves  measured  directly  over  the  column  of  interest.   All  sensors 
were  placed  within  a  2  meter  cube  of  water.   Simultaneous  time  series 
measurements  were  made  for  twenty-minute  durations.   The  experiments 
were  conducted  at  four  depths.   In  this  paper,  Run  Five,  at  a  depth  of 
7.3  meters,  will  be  analyzed  in  detail.   The  spatial  configuration  of 
the  sensors  is  shown  in  Figure  3.   The  actual  geometrical  separations 
are  shown  in  Figure  4. 

During  the  24  hour  period  of  the  experiment,  a  bi-modal  surface 
swell  of  about  0.3  meters  height  was  visually  observed  with  one  component 
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Carriage  System  Used 
to  Mount  Instruments 


Research  Tower  Schematic 
Figure  2. 
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having  a  period  of  about  7  seconds  from  285   true  and  the  other  compo- 
nent with  a  period  of  about  5  seconds  from  260   true.   The  surface  was 
observed  to  be  non-breaking  with  a  wind  capillary  system  present.   The 
wind  was  from  310   true  at  a  speed  of  about  8  knots.   A  moderate  temper- 
ature gradient  was  present  with  a  temperature  of  17.5  C  at  the  surface 
decreasing  to  15  C  at  a  depth  of  18  meters. 
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III.  INSTRUMENTATION 

A  brief  description  of  the  instruments  used  will  be  given  in  this 
section.   A  more  detailed  description  can  be  found  in  the  papers  of 
Bordy,  Smith,  Duchock,  Rautmann,  and  Seymour.   Instrument  details  are 
given  in  manufacturers'  technical  manuals. 

A.  SOUND  VELOCITY  SENSOR 

Sound  velocity  was  measured  with  a  Ramsay  MK-1  deep  sea  probe  util- 
izing only  the  sound  velocimeter  portion  of  the  probe.   The  Ramsay  veloc- 
imeter  is  of  the  "sing-around"  type  consisting  of  a  transducer  and 
transmitter  which  transmits  a  pulse  of  4  MHz  into  a  25  centimeter  sound 
path.   The  pulse  is  reflected  twice  to  reduce  errors  due  to  water  motion, 

B.  TEMPERATURE   SENSOR 

Temperature  was  measured  using  a  specially  constructed  Wheatstone 
bridge  circuit  utilizing  a  temperature ' sensitive  thermistor  as  one  leg 
of  the  circuit.   For  Run  Five,  two  such  circuits  were  employed.   The 
type  of  temperature  sensor  used  was  a  type  K  496  isocurve  thermistor 
manufactured  by  Fenwall  Electronics,  Inc.,  Farmington,  Massachusetts. 
The  specific  thermistors  were  extracted  from  an  expendable  bathythermo- 
graph.  It  is  a  double-bead  thermistor  which  incorporates  a  special 
aging  process,  giving  stable  and  repeatable  accuracies  not  usually  ob- 
tained with  ordinary  thermistors.   These  thermistors  have  a  time  con- 
stant of  150  milliseconds. 

C.  BAYLOR  WAVE  GAUGE 

The  Baylor  Model  13528  is  a  wave  profile  recording  system  which 
measures  instantaneous  sea  surface  elevation  as  a  function  of  time. 
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Measurements  are  accurate  within  one  percent  of  actual  height.   The  sys- 
tem consists  of  three  components:  the  wave  staff,  transducer  and  recorder, 
The  transducer  unit  measures  the  electrical  resistance  of  the  length  of 
the  staff  above  the  water  surface.   The  conducting  seawater  provides  a 
short  circuit  across  the  cables  (wave  staff)  immersed  in  the  water.   A 
direct  current  voltage  proportional  to  the  height  of  the  sea  surface  is 
taken  across  the  resulting  fluctuating  resistance. 

D.  WATER  CURRENT  METER 

The  water  particle  velocities  were  measured  with  an  Engineering 
Physics  Company  water  current  meter  Model  EMCM-3B.   The  EMCM-3B  is  an 
electromagnetic  current  measuring  device  that  measures  two  orthogonal 
components  of  water  particle  velocity  from  zero  to  five  meters  per  sec- 
ond.  Two  d.c.  voltage  signals  proportional  to  the  magnitude  of  the 
orthogonal  velocity  components  are  produced.   The  measured  response  is 
flat  to  at  least  0.5  seconds,  the  high  frequency  limit  of  the  calibra- 
tion facility.   The  maximum  output  error  is  one  percent  of  full  scale 
reading. 

E.  SALINITY  SENSOR 

The  salinity  sensor  employed  was  part  of  the  Bissett-Berman  STD 
Model  9006.   An  inductively-coupled  sensor  detects  conductivity,  which 
is  also  compensated  for  temperature  and  pressure  effects,  thereby  pro- 
ducing an  output  totally  dependent  upon  salinity.   The  time  constant  for 
the  salinity  sensor  is  .35  seconds. 

F .  MODULATION 

In  order  to  measure  environmentally  caused  amplitude  modulation  in 
the  ocean,  a  high  frequency  acoustic  carrier  wave  was  directed  at  a 

18 


receiving  hydrophone  mounted  along  the  acoustic  axis  at  a  separation  of 
two  meters.   An  envelope  detection  system  was  used  to  measure  the  ampli- 
tude modulation  of  the  received  signal.   A  USRD  type  F-27  transducer  and 
an  Atlantic  Research  Corporation  type  LC-10  hydrophone  were  used  in  the 
experiments . 

Go   SOUND  VELOCITY  (PHASE) 

In  determining  sound  velocity  by  phase  measurements,  one  F-27  trans- 
ducer was  used  along  with  one  LC-10  and  one  LC-5  hydrophone.   See  Rautmann 
L7]  for  further  details. 
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IV.   DATA  PRE-PROCESSING 

Run  Five  was  started  at  1728  and  ended  at  1748  on  21  October  1972. 
During  the  run,  data  was  recorded  on  a  Sangamo  Model  3500,  14  channel 
FM  tape  recorder.   The  data  consisted  of  temperatures  1  and  2,  salinity, 
wave  height,  u  and  w  water  particle  velocity  components,  sound  modula- 
tion, and  sound  velocity.   The  taped  data  was  transcribed  to  a  recti- 
linear eight  channel  Clevite  Brush  Recorder  after  amplification  by 
Krohn-Hite  Model  3322  variable  filters. 

A  Calraa  Company  Model  480  mechanical  digitizer  was  used  to  digitize 
the  continuous  strip  chart  information  at  100  points  to  the  inch.   This 
corresponds  to  a  sampling  interval  of  0.2032  seconds.   This  information 
was  recorded  onto  a  7  track  magnetic  tape  which  is  compatible  with  the 
CDC  6500  computer  at  the  Fleet  Numerical  Weather  Central,  Monterey.   The 
data  was  transferred  to  cards  using  program  CONVERT  (Lynch,  1970)  so  that 
it  could  be  read  for  analysis  by  the  IBM  360  at  the  Naval  Postgraduate 
School  which  uses  9  track  magnetic  tapes.  The  calibration  factors  used 
to  convert  the  data  into  the  appropriate  units  are  listed  in  Table  I. 
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TABLE  I 


CALIBRATION  FACTORS  FOR  DATA 

CALIBRATION  FACTOR 

MEASUREMENT 

UNITS 

.23 

Temperature 

°C 

.508 

Salinity 

ppt 

.58 

u,  w 

meters/second 

.34 

Sound  Velocity 

meters/second 

.211 

Fractional  Modulation 

.775 

Wave  Height 

meters 

This  information  is  included  here  to  facilitate  future  analysis  of  this 
data. 
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V.   ANALYSIS  OF  DATA 

A.   FOURIER  FILTER 

The  experiment  was  designed  to  study  processes  in  the  approximate 
frequency  band  0.02  to  2  Hz  corresponding  to  time  scales  of  0.5  seconds 
to  50  seconds.   The  initial  analysis  of  the  raw  data  showed  that  the 
energy  density  spectra  of  some  parameters  had  excessive  leakage  from 
the  very  low  frequencies.   The  time  correlations  were  of  the  order  of 
minutes.   Any  effects  caused  by  forcing  functions  outside  the  frequency 
band  of  interest,  such  as  tidal  or  internal  wave  effects,  will  not  be 
considered.   Thus  it  is  necessary  to  eliminate  the  low  frequency  com- 
ponents in  the  raw  data  records.   This  was  accomplished  by  using  a 
Fourier  Filter. 

The  first  step  is  to  generate  a  Fourier  series  of  the  record.   The 
fast  Fourier  Transform  algorithm  is  used  to  conserve  computer  time: 

00 

An    V' 

F(t)  =  ~  +  /  (A  cos  niut  +  B  sin  nit) 
2    L>       n  n 

n=l 
where       F  =  parameter 

t  =  time 

njj  =  nth  harmonic  of  fundamental  radial  frequency  uo 

Next,  a  time  series  is  generated  using  only  the  low  frequency  components 
that  are  to  be  eliminated.   This  time  series,  composed  only  of  the  low 
frequency  components,  is  then  subtracted  from  the  original  time  series 
for  each  value,  and  the  remainder  is  the  filtered  value  of  the  parameter 
at  each  given  time.   In  this  experiment,  the  first  fifteen  low  frequency 
components  were  used  to  generate  the  filtered  time  record,  corresponding 
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to  the  components  with  a  fifty-five  second  periodicity  and  greater.   Al- 
though certain  records  needed  filtering  while  others  did  not,  all  records 
were  filtered  for  consistency. 

The  effect  of  the  Fourier  Filter  on  the  temperature  2  spectra  can  be 
seen  in  Figure  5.   Notice  the  removal  of  the  energy  from  0  to  .018  Hz  in 
the  energy  density  spectrum.   After  this  point,  the  energy  spectra  coin- 
cide for  the  remainder  of  the  spectrum.   Also  of  importance  is  that  the 
phase  and  coherence  spectra  remain  unaltered  in  the  region  of  interest, 
greater  than  .018  Hz. 

B.   SPECTRAL  ANALYSIS 

The  energy-density  spectrum  is  calculated  from  the  auto-covariance 
function  of  the  recorded  time  series: 

.   t+T 

flu  (t)  =  lim  -J   x  (t)x  (t  +  t)  dt 

T-oo     t 

where  the  auto-covariance  function  is  expressed  as  a  function  of  the  lag 
time,  t.   If  the  auto-covariance  is  normalized  by  division  by  the  vari- 
ance, it  is  called  the  auto-correlation  function.   The  Fourier  transform 
is  the  power  or  energy-density  spectrum. 

T/2 

s       ,<-\         ,.         1'"  u       /x       -i2rrfT     , 

*lx  (f)   =   lim    -  J  0ii  (t)    e  dT 

T-°°    -T/2 

A  Parzen  lag  window  was  applied  to  the  functions  to  account  for  the 
finite  length  of  the  records  in  computing  the  spectra.   Since  the  Parzen 
window  has  no  negative  side  lobes,  this  eliminates  any  numerical  insta- 
bility problems  that  could  arise  in  the  analysis. 

Let  P(t)  denote  the  Parzen  Window  Function,  then  the  modified  energy 
density  spectrum  can  be  written: 
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-i2TTfT 


$11  (f)    =    J       P(T)01l(T)e"1Z    tT    dT 


The  energy-density  spectrum  is  the  statistical  estimate  of  the  energy- 
density  distributed  over  the  frequency  bands  of  the  spectrum.   The  total 
potential  energy-density  of  the  sea  surface  height  is  proportional  to  the 
total  area  under  the  wave  height  spectrum.   Similarly,  the  area  under 
the  velocity  spectrum  is  proportional  to  the  kinetic  energy-density.   The 
area  under  each  curve  is  equal  to  twice  the  variance  of  the  random  vari- 
able. 

The  cross-spectral  density  is  defined  as  the  Fourier  transform  of 

the  cross-covariance  function.   The  cross-covariance  function  is  similar 

to  the  auto-covariance  function  except  that  a  second  time  series  record 

is  lagged  against  the  first  instead  of  a  single  record  lagged  against 

itself.   The  cross-covariance  function  indicates  the  dependence  the 

values  of  one  record  have  on  that  of  another  for  various  lag  times. 

The  cross-covariance  function  is: 

1  J/2 
0is(T)  =  lim  -J    x  (t)x  (t  +  t)  dt 
T-°°    -T/2 

The  modified  cross  spectral  density  is: 

00 

§12  (f)  =  J  POO0i2(t)  e"i2nfT  dT 

-co 

The  cross-spectrum  can  be  defined  in  terms  of  its  real  and  imaginary 
parts : 

§12  (f)  =  Ci2(f)  -  iQia(f) 
where  the  real  part  of  §12  (f)  is  known  as  the  co-spectrum: 

0D 

Cis(f)    =    2 J    [§12  CO    +    il3(-T)]     COS     (2TTfT)     dT 
0 
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and  the  imaginary  part  of  $ig  (f)  is  the  quadrature  spectrum: 

00 

Qia(f)  =  2j'  [$12(T)  +  ?i2(-t)]  sin  (2nfT)  dT 
0 

The  spectral  phase  angle  represents  the  average  angular  difference 

by  which  the  cross  correlated  components  of  x9  (t)  lead    those  of  x..  (t) 

in  each  spectral  frequency  band.   The  phase  angle  is  defined  by: 

el2  (f )  =  arctan  g£^fi] 

e^s  (f)  ranges  from  -18CP  to  +180°  . 

The  coherence  function  is  defined  by: 

Yl2(f)  =   tSir(f)]2 
Yl2U;     *li(f)*22(f) 

with  the  provision  that: 

I  *i2(f)|  *  *n(f)*22(f) 

which  results  in: 

0  *  Y12  (f)  ^  1 

The  coherency  spectrum  is  a  measure  of  the  correlation  between  two  time 
series  as  a  function  of  frequency.   Three  possible  causes  of  coherency 
values  being  less  than  unity  are:   the  system  is  not  linear,  excess  noise 
is  present  in  the  measurement,  or  the  output  time  series  is  due   to  the 
input  time  series  as  well  as  other  functions  (Bendat  and  Piersol  [_l]). 

C.   SPATIAL  CORRELATION 

The  data  were  measured  as  a  time  series  record  at  specific  points  in 
space.   The  majority  of  theory  in  this  area  of  study  requires  knowledge 
of  spatial  correlations.   Hence,  meaningful  comparisons  cannot  be  made 
without  measures  of  spatial  correlations.   The  operation  described  below 
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is  designed  to  convert  the  time  domain  record  to  a  spatial  one.   The  tem- 
poral measurements  are  transformed  to  a  Lagrangian  frame  knowing  the  ve- 
locity and  acceleration  at  the  point  of  parameter  measurement. 

It  is  assumed  a   priori  that  the  correlation  lengths  are  small  com- 
pared to  the  spatial  scales  of  motion.   At  7.3  meters  depth  the  orbital 
excursion  of  a  10  second  wave  having  a  rms  wave  height  of  30  centimeters 
(mean  conditions  during  experiment)  is  80  centimeters.   It  is  further 
assumed  that  the  temperature  field  is  isotropic  to  simplify  the  analysis. 
Measurements  of  velocity  in  the  X  and  Z  plane  were  taken  simultaneously 
along  with  temperature.   Although  the  velocity  components  were  measured 
a  small  distance  from  the  thermistor  it  is  consistent  to  assume  that  the 
velocity  components  can  be  transferred  to  the  thermistor  location,  r  , 
due  to  the  small  separation  distance  involved.   If  one  now  considers  the 
thermistor  moving  with  the  velocity,  U,  and  acceleration  components,  it 
will  travel  a  path  in  the  X-Z  plane.   At  each  instant  the  velocity  is 
measured,  the  distance  traveled  by  the  thermistor,  r,  is  calculated  by 

2 

r  =  r  ^ +  u&t  (C-l) 

o    2At 

(At=  time  increment) 
whereupon  the  next  temperature  is  assigned  to  this  distance  "r."   Thus 
an  array  is  created  where  a  temperature  is  assigned  for  a  certain  dis- 
tance "r"  in  the  X-Z  plane.   For  each  temperature  measurement  in  time 
there  is  an  associated  distance  giving  T(r,t).   In  order  to  facilitate 
calculation  of  the  auto-correlation  function,  any  holes  in  the  array 
were  filled  by  linearly  averaging  the  temperature  over  the  distance 
increment . 

The  question  arises  as  to  the  validity  of  the  concept  in  relation  to 
the  occurence  of  a  negative  velocity  which  would  necessitate  a  backward 
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movement  of  the  thermistor,  bringing  about  the  possibility  of  two  differ- 
ent temperatures  being  assigned  to  the  same  point  in  the  X-Z  plane.   The 
probability  of  this  happening  is  small,  and  the  effect  of  its  occurrence 
is  negated  by  averaging  the  first  temperature  with  the  second  one.   If 
by  chance  a  third  temperature  were  assigned  this  point,  it  would  be  aver- 
aged with  the  previously  averaged  temperature;  hence  the  last  temperature 
would  always  have  a  weighting  factor  of  50  per  cent. 

A  similar  scheme  using  an  Eulerian  frame  can  be  derived  by  expanding 
the  temporal  measurement  at  a  point  in  a  Taylor  series: 

m/    x    „./    n    rAT(r0,t)   ,   r^A^ T(r0  ,  t)   ,  ,„  oN 

T(r,t)  =  T(rQ,t)  +  ^  '  +     2Ar~°       "•        (C'2) 

where  Ar  =  r  -  r  and  Ar  is  small.   Employing  Taylor's  hypothesis  that 

Ar 
over  short  distances  and  time   U  =  —  ,  Equation  (C-2)  can  be  rewritten: 

Tfr  t)  -  Tfx  t)  -»■  ^^^  +  r^TCr^t) 

The  two  schemes  are  equivalent,  but  computationally  the  Lagrangian  form 
is  more  convenient  and  is  employed  here. 

D.   AUTO-CORRELATION  FUNCTION 

With  the  conversion  of  a  time  series  record  to  a  spatial  record  in 
the  preceding  section,  the  description  of  the  parameter  in  the  medium  by 
the  spatial  auto-correlation  function  is  now  possible.   The  importance  of 
the  spatial  distribution  of  the  parameter  is  illustrated  in  the  following 
example.   A  sharp  temperature  boundary  occuring  within  an  acoustical 
wavelength  of  the  propagated  sound  may  cause  reflection.   The  linear 
dimensions  of  the  temperature  inhomogegeity  may  cause  scatter,  reflection, 
or  focusing  of  sound  waves.   The  physical  features  of  a  parameter  may  be 
characterized  by  the  auto-correlation  function. 
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The  computation  of  the  spatial  auto-correlation  function,  0(p,t), 
from  the  temperature  record,  T(r,t)  is  defined: 

,   L/2 
0(p3t)=Tr   T(r,t)'T(r+p,t)  dr  (D-l) 

L  -L/2 

0(0, t) 

where  p  is  the  lag  distance  and  L,  the  total  length  of  the  field.   Fig- 
ure 6  shows  the  spatial  auto-correlation  function  for  temperature  1  and 
2.   Figure  7  shows  the  time  auto-correlation  function  for  the  same  para- 
meters.  Note  the  rapid  decay  to  zero  by  both  functions  and  also  the 
sinusoidal  components;   thus  both  cases  reveal  the  presence  of  periodic 
components,  namely  wave  induced  particle  motions. 

The  distribution  of  temperature  2  in  time  is  shown  to  closely  ap- 
proximate a  Gaussian  distribution  in  Figure  8,  Figure  9  (using  the 
Fourier  Filter).   Thus  a  good  description  for  the  function  over  short 
time  lags  would  be: 

e 

0(t)  =  e~  to2  (D-2) 

where  t   is  the  correlation  time.    Assuming  a  similar  spatial  Gaussian 
distribution  gives: 

0(P)  =  e  po2  (D-3) 

The  mean  "size"  of  the  inhomogeneities  might  reasonably  be  described 
by  the  value  of  p  at  which  the  temperature  0(p)  drops  to  1/e  denoted  by 
po  in  Equation  D-3.   These  distances  are  known  as  correlation  distances 
(Liebermann  [4]).   For  temperature  records  1  and  2,  their  values  of  p0 
were  26  and  27  cm  respectively. 

The  consistent  results  for  correlation  distances  show  that  the  con- 
version from  time  to  the  spatial  domain  is  reasonable.   Hence,  the 
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FIGURES.     TEMPERATURE    2    PROBABILITY    DENSITY    FUNCTION 
WITHOUT    FOURIER    FILTER 
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Figure  9.   Temperature  2  Probability  Density  Function  With  Fourier  Filter. 
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approximate  mean  size  of  the  temperature  inhomogeneities  at  a  depth  of 
7.3  meters  are  characterized  at  about  26  centimeters.   The  difference  in 
correlation  times  for  temperature  1  and  2  as  seen  in  Figure  7  (4.5  sec- 
onds versus  1.4  seconds)  could  be  due  to  higher  turbulence  in  the  loca- 
tion of  thermistor  2  which  would  result  in  smaller  inhomogeneities.   An 
alternative  possibility  would  be  a  malfunction  of  the  thermistor.   Regard- 
less, the  conversion  of  a  time  series  to  a  spatial  domain  has  the  effect 
of  averaging  out  the  localized  turbulence  by  filling  the  spatial  tempera- 
ture array  with  smoothed  values  of  temperature  variations.   Thus  the 
spatial  processing  tends  to  average  by  minimizing  the  effect  of  localized 
turbulence . 

A  satisfactory  distribution  function  for  the  near  surface  zone  must 
take  into  account  periodic  motion  induced  by  surface  waves.  Therefore, 
a  more  accurate  description  of  the  auto-correlation  function  would  be: 


0(P,t)  = 


P2 


Po   i  ■?-  +  £i[A  cos  n(kx-at)  +  B  sin  n(kx-at)] 
JL2   n=  1  n  n  J 


However,  this  function  is  entirely  too  cumbersome.   R(p)  =  e   p0   will  be 
used  because  it  describes  the  function  very  well  through  the  correlation 
distance  up  to  the  first  zero  crossing  and  is  capable  of  being  integrated 
which  is  important  in  a  later  section. 

E.   WILSON'S  EQUATION 

Wilson's  Equation  [2]  was  originally  developed  to  relate  the  speed 
of  sound  to  the  properties  in  ocean  water  that  affect  it,  namely  pressure, 
temperature,  and  salinity.   An  abbreviated  form  of  Wilson's  Equation  is: 

c  =  1449  +  4.62T  -  0.055T2  +  0.0003T3  +  (1.39- . 012T)  (s-35)  +  0.017z 
where   c  -    speed  of  sound  in  meters  per  second 
T  -  temperature  in  degrees  Centigrade 
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s  =  salinity  in  parts  per  thousand 
z  -  depth  in  meters 

One  of  the  initial  aims  of  the  experiment  was  to  show  that  the  vari- 
ance of  the  speed  of  sound  was  equal  to  the  sum  of  the  variances  of  its 
constituents  multiplied  by  their  respective  coefficients.   This  would 
establish  the  validity  of  Wilson's  Equation  in  terms  of  variances  as 
well  as  indicate  the  accuracy  of  the  temperature,  salinity,  and  depth 
data. 

Wilson's  Equation  applied  in  this  manner  was  not  successful  for  the 
experiment.   Let 

var  c  =  (  (c-c)/ 
dT  -  <  T-T  > 
ds  =  (  s-s  ) 
where  the  overbar  implies  averaging  and  (  )  indicates  an  ensemble  average, 
After  simplifying  by  retaining  only  first  order  terms  of  temperature  and 
salinity  and  dropping  the  small  depth  term, 

var  c  -  <(5.04dT+  1.39  ds)2> 
var  c  =  <(25.40dT2  +  14.00  dT»ds  +  1.93ds2)> 
and  dT'ds  =  covariance  (T*s)  =  cov  (T«s). 

In  order  to  measure  the  variance  of  c  using  Wilson's  Equation,  both 
temperature  and  salinity  must  be  measured  at  the  same  point  at  the  same 
time.   The  time  constant  for  the  salinity  sensor  was  .35  seconds,  and 
for  the  thermistors  .15  seconds.   The  average  spatial  distance  between 
the  instruments  was  2.14  meters. 

Correlation  distance  has  been  described  in  a  previous  section  as 
representative  of  the  mean  size  of  the  inhomogeneiety .   In  order  to  en- 
sure that  temperature  and  salinity  sensors  provide  accurate  results, 
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they  must  be  located  within  the  same  inhomogeneity  along  with  the  sound 
velocity  meter.   This  is  obtained  only  if  all  the  instruments  are  placed 
within  the  smallest  correlation  distance  of  the  parameters  involved.   As- 
suming that  the  mean  size  of  a  salinity  inhomogeneity  would  be  no  smaller 
than  that  for  temperature,  then  all  instruments  would  have  to  have  been 
placed  within  26  centimeters  in  this  experiment  (as  previously  determined 
from  temperature  correlation  distances)  in  order  to  provide  an  accurate 
evaluation  of  Wilson's  Equation.   For  the  experimental  data  taken,  an 
application  of  Wilson's  Equation  in  this  case  would  be  futile  as  Seymour 
[8]  showed  in  his  computer  results. 

As  the  distance  between  instruments  increases  in  the  random  medium, 
one  would  expect  that  local  mechanisms  operating  within  one  correlation 
distance  to  bear  less  and  less  resemblance  to  mechanisms  operating  at 
increased  correlation  distances;  hence,  the  temperature  and  salinity 
would  approach  statistical  independence  with  increasing  distance,  causing 
cov  (T's)  to  approach  zero  by  definition.   This  was  found  to  be  the 
case  in  Run  Five  where  the  cov  (T«s)  integrated  over  the  spectrum  was 
found  to  be  -.8  x  10   (C  -PPT).   The  variance  of  temperature  was  1.6  x 

—  A      r\     0  0  0 

10    (C  )  ,  the  variance  of  salinity  0.22  x  10    (PPT)  ,  and  the  variance 

-4  2 

of  sound  velocity  1.1  x  10    (meters/second)  . 

The  importance  of  the  waves  to  the  variance  of  the  parameters  is 

indicated  by  band-pass  filtering.   Band-pass  filtering  was  accomplished 

by  high-pass  filtering  with  the  Fourier  filter  described  before  and 

calculating  the  variance  by  summing  the  spectral  contributions  only  of 

frequencies  up  to  0.16  Hz.   The  band-pass  filter  was  chosen  because  it 

included  that  part  of  the  variance  where  organized  wave-induced  particle 

motion  was  dominant  in  its  contribution  to  the  variance,  and  the  effects 


36 


of  turbulence  and  instrument  noise  were  insignificant.   Table  II  shows  the 
results  of  the  band-pass  filter. 

TABLE  II 

VARIANCES  OF  SOUND  VELOCITY,  TEMPERATURE ,  AND  SALINITY 

Variance  Entire  Spectrum      Band-pass  Filter     %  Variance 

in  wave  dom- 
inated region 

var  c  1.1  x  10*4  1.0  x  10~4  91% 

var  T  1.6  x  10~4  0.6  x  10~4  38% 

var  s  21.5  x  10~4         19.1  x  10~4  90% 

On  inspection,  the  high  value  of  variance  for  salinity  makes  it 
questionable  as  one  would  expect  a  value  of  the  same  order  of  magnitude 
as  the  temperature,  not  13  times  as  great.   Further  investigation  in- 
dicates that  the  measured  conductivity  is  compensated  continuously  for 
pressure  and  temperature  changes  with  the  result  that  the  indicated  sal- 
inity has  a  tendency  to  oscillate;  thus  the  high  value  for  the  variance 
of  salinity  should  be  disregarded.   This  is  significant  in  itself,  and 
when  combined  with  a  lack,  of  the  covariance  term,  cov  (T*s),  it  makes 
any  proof  of  Wilson's  Equation  only  conjecture.   Hence,  future  investi- 
gations must  use  more  appropriate   salinity  instrumentation  in  the  near 
surface  region  of  the  ocean. 

The  frequency  of  the  sound  velocimeter  was  4  MHz.   In  this  frequency 
range,  the  bubbles  present  in  the  near  surface  zone  have  no  effect  on  sound 
propagation  due  to  the  extremely  small  wavelength  of  the  sound  wave. 
Thus  by  measuring  the  variance  of  the  sound  velocimeter  and  comparing  it 
by  Wilson's  Equation  to  the  variances  of  the  temperature, salinity  and 
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their  respective  covariance,  one  gets  a  crosscheck  on  the  refraction  of 
sound  in  the  medium  due  only  to  salinity  and  temperature. 

F.   AMPLITUDE  MODULATION  IN  A  RANDOM  MEDIUM 

In  an  ideal  homogeneous  undisturbed  medium,  a  sound  wave  propagated 
at  a  fixed  frequency  and  amplitude  arrives  in  the  far  field  undiminished 
except  for  spherical  spreading  and  absorption  losses.   An  adequate  de- 
scription of  the  process  for  a  homogeneous  medium  is  given  by  the  general 

wave  equation: 

2 
V2P  =  1/c2  ^~  (F-l) 

Bt 

2 
where       V  =  LaPlacian  operator 

P  -  pressure 

c  -  speed  of  propagation  in  medium. 
The  medium  under  study  is  the  near  surface  ocean.   Bubbles,  tempera- 
ture and  salinity  patches,  and  scatterers  are  all  present  in  a  random 
distribution.   Compounding  the  problem  is  the  fact  that  the  inhomogene- 
ities  are  in  constant  motion  due  to  wave  particle  motion  or  convective 
forces.   As  the  sound  wave  moves  through  a  medium,  the  inhomogeneities 
cause  a  variation  in  the  speed  of  sound  which  results  in  refraction.   If 
the  inhomogeneities  were  motionless,  reflection  and  refraction  of  a 
propagating  sound  beam  would  cause  the  amplitude  of  the  sound  beam  to 
be  constant  at  a  given  receiver.   If  the  inhomogeneities  have  freedom 
to  move  in  and  out  of  the  path  of  the  sound  beam,  then  the  amplitude  of 
the  sound  beam  at  the  receiver  will  fluctuate.   Hence,  the  wave  equa- 
tion must  be  modified  to  include  spatial  and  temporal  changes  in  the 
speed  of  sound.   The  change  in  the  index  of  refraction  with  time  and 
space  is  defined: 
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(-  t)   .  SlL^l  =   Co  -f-  dc(r,t)   m    x  +  dc(r,t)   =    ±  +  a 

C<  Cq  Cq 

,  „  _   dc(r.t) 

where:  <x  = 

co 

u.   -   Index   of    refraction 

Cq=  Average   speed   of   sound 

c  =   co+dc~  Local   speed   of    sound 

ck  =  RMS   refractive   index  variation   from  unity 

t  -  Time 

r  -   Position  vector 

Previous  measurements  of  Oi   for  temperature    inhomogcncities  have 

-A  r     I 

yielded  values  such  as  1.18  x  10    (Stone  and  Mintzer  in  Laboratory  LlOj) 
and  7.0  x  10  "  (Liebermann  L6]  in  Ocean  Experiments  at  depths  of  50 
meters).   In  Run  Five  of  this  experiment,  the  calculated  a   for  tempera- 
ture inhomogeneities  at  a  depth  of  7.3  meters  was  3.0  x  10   .   This  was 
calculated  by  taking  the  first  derivative  of  Wilson's  Equation  with  re- 
spect  to  temperature  and  substituting  the  rms  value  of  the  temperature 
variance  into  only  temperature  terms  up  to  the  third  order;  the  average 
speed  of  sound  was  calculated  using  Wilson's  Equation  with  time  averaged 
values  for  temperature,  salinity,  and  pressure.   In  comparing  Lieber- 
mann 's  results  to  Run  Five,  the  implication  is  that  the  temperature  in- 
homogeneities are  smaller  near  the  surface.   This  seems  reasonable  as 
water  particle  motion  caused  by  wave  action  decays  exponentially  with 
depth;  hence  better  mixing  would  be  expected  nearer  the  surface  reducing 
the  size  of  the  temperature  inhomogeneities. 

In  their  laboratory  studies  of  amplitude  modulation  in  a  random  med- 
ium, Mintzer  and  Stone,  [6j,  defined  the  coefficient  of  variation,  cv: 

2  =  ,«P2>  -<P>-)  (F_3) 


"<py 
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where,      cv  -  The  coefficient  of  variation 

P  =  The  amplitude  of  the  received  signal 

The  coefficient  of  variation  as  defined  by  Mintzer  and  Stone  is  related 
to  the  fractional  modulation  measured  in  Run  Five  by  the  following  equa- 
tion : 

2 
(.0253) (cv)   =  var  Modulation 

The  development  of  this  relationship  is  found  in  Appendix  I  of  Smith  1.9]. 

The  majority  of  theory  dealing  with  propagation  in  a  moderately  in- 
homogeneous  medium  employs  the  single  scattering  approximation  to  the 
wave  equation.   Mintzer  and  Stone  l6]  have  shown  that: 

CO 

(cv)2   =   <*2Ko2b  J  0(p)    dp  (F-4) 

o 
where: 

K     -  Wave   number 
o 

b  -    range  from  source  to  receiver 

0 (p )   =  the  correlation  function  taken  along  the  path 
from  source  to  receiver. 

ex     =   root  mean  square  index  of  refraction  variation  from 
unity. 

Previously,  it  was  determined  that  the  temperature  auto-correlation 
was  described  by  0 (p )  =  e   p0  .   Other  parameters  affecting  the  index 
of  refraction  showed  similar  Gaussian  distributions;   therefore,  sub- 
stituting the  expression  for  0 (p )  (representing  a  combined  temperature, 
salinity  and  bubble  inhomogeneity  auto  correlation  function)  into  equatiion 
(F-4)  and  integrating,  one  gets  the  following  result: 

(cv)2  =  %  /ii   K  2a2p  r  (F-5) 

o    o 

In  Run  Five  the  sound  modulation  was  obtained  for  a  sound  wave  of 

frequency  60  kHz.   This  frequency  has  a  wavelength  of  2.5  centimeters 
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which  happens  to  correspond  to  a  peak  in  the  resonant  bubble  regime   of 
the  near  surface  zone  at  a  depth  of  7.3  meters  in  this  experiment,  as 
determined  from  Rautmann  [7].   By  using  two  hydrophones  to  determine 
phase  change  over  a  78  cm  path,  Rautmann  converted  the  phase  change  to 
a  change  in  the  speed  of  sound.   At  7.3  meters  depth  during  the  same 
18  hour  time  period,  an  interpolated  value  for  60  kHz  from  Rautmann 's 
tables  indicates   a  c(measured)  -  c (theoretical)  =  -1  m/sec.  C(measured) 
is  the  speed  of  sound  based  on  phase  measurements  while  c (theoretical)  is 
taken  from  Wilson's  Equation  using  average  temperature,  salinity,  and 

depth.   Thus  the  variance  of  the  speed  of  sound  at  60  kHz  was  determined 

2 
to  be  1  (m/sec)   due  to  temperature,  salinity,  and  bubble  inhomogeneities . 

Converting  this  to  an  ot   value  gives  the  relatively  high  value  of 

-3 
.67  x  10   .   This  value  is  based  on  bubbles,  temperature,  and  salinity, 

not  just  temperature  alone  which  gave  a  value  of  .25  x  10   .   Since  the 
salinity  values  were  highly  suspect  and  along  with  the  lack  of  a  cov(T*s) 
term,  one  cannot  accurately  determine  the  exact  effect  of  bubbles  on  or. 
However  if  one  makes  the  modest  assumption  that  the  salinity  contribution  to 
dT'dswill  be  equal  to  the  temperature  contribution  or  less,  then  the 
temperature  and  salinity  contribute  less  than  one  percent  of  the  value 
of  <y .   This  shows  the  enormous  importance  of  bubbles  in  the  near  surface 
region  upon  the  refraction  of  sound  at  60  kHz. 

Further  proof  of  the  importance  of  bubbles  can  be  seen  in  the  spec- 
tral analysis  of  wave  height  versus  temperature,  salinity,  and  sound  mod- 
ulation (see  Figures  10,11,  and  12).   In  the  coherence  and  phase  spectra 
of  wave  height  versus  temperature  and  salinity,  small  consistent  values 
of  coherence  and  phase  are  found  in  the  wave  dominant  portion  of  the 
spectrum  (.05  to  .22  Hz).   However,  for  wave  height  versus  sound  modulation, 
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the  coherence  value  approaches  unity  and  the  phase  is  constant  in  the 
wave  dominant  region.   This  implies  that  wave  height  is  linearly  re- 
lated to  sound  modulation.   Since  the  effect  of  the  change  in  pressure 
due  to  wave  height  is  insignificant  on  the  sound  modulation  at  this  depth 
then  the  fluctuation  in  modulation  must  be  due  to  the  variations  in  the 
inhomogeneities.   With  the  coherence  of  wave  height  and  temperature  and 
wave  height  and  salinity  being  low,  the  process  of  elimination  leaves 
only  one  type  of  inhomogeneity  --  bubbles.   This  supports  the  contention 
for  this  experiment  that  the  dominating  cause  of  change  in  the  speed  of 
sound  was  due  to  bubbles  as  evidenced  by  very  high  coherence  of  wave 
height  and  modulation. 

This  high  coherence  can  be  due  to  two  wave  related  mechanisms:   1)  the 
bubbles  (inhomogeneities)  are  moved  in  and  out  of  the  sound  path  by  the 
wave  induced  motion   and,  2)  the  bubbles  are  "squeezed"  by  the  wave  induced 
pressure  changes  resulting  in  a  change  in  their  resonance  characteristics. 

In  applying  the  data  from  Run  Five  to  Mintzer's  equation,  "Po>"  the 
mean  radius  of  the  scattering  inhomogeneity,  is  found  to  be  equal  to  8.8 
centimeters.   This  is  not  to  say  that  an  inhomogeneity  exists  with  a 
radius  of  8.8  centimeters,  but  the  effect  of  bubbles,  temperature,  and 
salinity  provides  an  "effective"  radius  of  8.8  centimeters  for  the  in- 
homogeneity. 

G.   THE  EFFECTS  OF  SURFACE  WAVES  ON  THE  MEASURED  PARAMETERS 

There  are  three  contributors  to  the  energy  density  spectrum  of  hori- 
zontal or  vertical  particle  motion:   wave  induced  particle  motion,   mean 
convective  motion,  and  turbulence.   Mean  convective  motion  was  essentially 
zero  for  the  particle  velocity  in  the  vertical  plane.   This  leaves  wave 
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induced  particle  motion  and  turbulence  as  the  two  dominant  sources  in 
supplying  energy  to  the  power  spectrum  of  vertical  particle  motion. 

Turbulence  is  random  three-dimensional  motion  which  can  dissipate 
mechanical  energy  to  internal  energy  through  a  cascade  of  eddies  of 
diminishing  size.   The  process  of  turbulent  cascade,  resulting  in  dis- 
sipation, includes  nonlinearity  and  three  dimensionality  as  well  as 
rotationali ty .   In  contrast  to  this  is  wave  induced  particle  motion 
which,  while  being  random  and  three  dimensional,  is  approximately  ir- 
rotational  and  nondissipative  (Lumley  and  Panofsky  [5]).   Figure  13 
clearly  establishes  the  irrotationality  of  the  wave  induced  velocities 
by  showing  an  almost  constant  phase  difference  of  90  between  u  and  w 
components  in  the  region  of  significant  energy  density.   The  90  phase 
difference  is  predicted  by  solving  the  linear  first  order  wave  equation 
in  terms  of  velocity  potential.   A  qualifying  condition  for  a  velocity 
potential  is  the  existence  of  irrotationality  in  the  fluid. 

The  horizontal  component  of  water  particle  velocity,  u,  will  not  be 
considered  in  terms  of  absolute  or  quantitative  values  as  the  velocity 
meter  measured  only  one  component  of  the  horizontal  velocity  during 
Run  Five. 

The  measured  energy-density  spectrum  of  the  vertical  water  particle 
velocity  has  been  plotted  in  Figure  14.   Figure  14  also  shows  the  cal- 
culated values  of  the  vertical  velocity  energy-density  spectrum  predicted 
using  linear  wave  theory  (Bordy  [2]): 

,2 


I  ,ni2   T  f   sinh  k(h+z)  1  \vy/e\ 
lw(f)l   =  L  2^    sinh  kh    J  l^(f) 


The  term  in  brackets  is  the  first  order  vertical  velocity  amplitude  - 
response  function  squared.  T\   is  a  measure  of  wave  height  as  a  function 
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of  frequency.  The  measured  velocity  was  equal  to  the  theoretical  value 
for  the  region  of  significant  energy-density.  Close  examination  of  the 
spectra  indicates  that  the  measured  and  theoretical  vertical  velocities 
begin  to  diverge  at  a  frequency  of  about  0.2  Hz.  This  marks  the  begin- 
ning of  the  transition  regime  to  turbulence.  At  0.3  Hz  turbulence  is  the 
dominant  process  in  the  energy  density  spectrum. 

Spectra  of  other  parameters  compared  to  waves  are  shown  in  Figures 
10,  11,  12,  and  15.   The  contribution  to  the  energy-density  spectra  of 
each  by  wave  induced  particle  motion  in  the  frequency  range  .05  Hz  to 
.22  Hz  approaches  95%  of  the  total  energy  with  the  exception  of  tempera- 
ture. 

The  phase   and  coherence  spectra  of  the  wave  height  versus  modula- 
tion, shows  the  coherence  approaching  unity.   This  strong  linear  rela- 
tionship of  waves  and  sound  modulation  was  attributed  to  bubbles  in  an 
earlier  section.   The  values  of  coherence  are  less  for  salinity,  sound 
velocity,  and  temperature,  with  the  least  coherence  for  wave  height  ver- 
sus temperature.   This  could  indicate  that  t lie  temperature  scalar  field 
is  more  affected  by  turbulence  and  thermal  dif f usivity*,  or  it  could  in- 
dicate an  instrumental  inadequacy.   Regardless,  it  is  noted  that  the 
temperature  has  a  significantly  less  linear  relationship  with  wave  height 
than  the  other  parameters  and  accordingly  only  about  407o  of  its  energy 
comes  in  the  wave  dominant  portion  of  the  spectrum. 


*  The  thermal  diffusivity  is  approximately  10  times  greater  than  salt 
diffusivity  in  the  ocean.   Thus  temperature  gradients  are  expected  to 
be  erased  much  quicker  than  salinity  gradients. 
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VI.   CONCLUSIONS 

A.  The  Fourier  high  pass  filter  proved  to  be  an  effective  method  of 
eliminating  low  frequency  components  from  the  signal  record  without  al- 
tering phase,  coherence,  or  energy  relationships  in  the  region  of  interest. 

B.  The  spatial  correlation  distances  for  temperature  determined  for 
measurements  at  a  point  were  established  to  be  about  25  centimeters. 

C.  Future  experiments  attempting  to  verify  Wilson's  equation  by  summing 
the  variances  of  its  constituents  (T,s,P)  should  require  the  instruments 
to  be  placed  within  a  correlation  distance;  this  insures  that  the  local 
mechanisms  affecting  temperature  and  salinity  inhomogeneities  are  iden- 
tical and  will  give  a  valid  covariance  term.   Further,  better  salinity 
instrumentation  is  required  for  the  time  scales  considered  in  the  near 
surface  region. 

D.  Sound  modulation  of  a  60  kHz  sound  beam  traversing  a  2  meter  path 
length  underwent   a  mean  attenuation  of  about  57,.   The  correlation  dis- 
tance of  an  effective  inhomogeneity  due  to  the  combined  temperature,  salin- 
ity, and  bubble  effects  was  found  to  be  8.8  centimeters  by  applying 
Mintzer  and  Stone's  Equation. 

E.  The  change  in  the  index  of  the  speed  of  sound  at  60  kHz  (a  resonant 
bubble  region)  was  found  to  be  997,  induced  by  bubbles. 

F.  The  energy  density  spectra  of  all  parameters  were  dominated  by  wave 
induced  motion  at  7.3  meters  depth. 
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G.   High  coherence  is  found  between  wave  induced  water  particle  motion 
and  sound  modulation  and  decreasing  coherence,  between  water  particle 
motion  and  salinity,  sound  velocity,  and  temperature. 

H.   First  order  wave  theory  accurately  predicts  orbital  water  particle 
velocities . 

I.   Wave  induced  water  particle  velocities  are  irrot ational. 
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ments are  analyzed  for  a  twenty  minute  record  taken  at  a  depth  of  7.3  m. 
The  spatial  correlation  distance  for  temperature  was  found  to  be  about  25 
cm.   Sound  modulation  of  a  60  kHz  sound  beam  underwent  a  mean  attenuation 
of  about  5%.   The  correlation  distance  of  an  effective  inhomogeneity  due 
to  the  combined  temperature,  salinity,  and  bubble  effects  was  found  to  be 
8.8  cm.   The  change  in  the  index  of  refraction  of  the  speed  of  sound  at 
60  kHz  (a  bubble  resonant  region)  was  found  to  be  99%  induced  by  bubbles. 
The  energy  density  spectra  were  dominated  by  wave  induced  motion.   First 
order  wave  theory  accurately  predicts  orbital  water  partical  motion  and 
its  irrotationality.   High  coherence  is  found  between  wave-induced  water 
particle  motion  and  sound  modulation  and  decreasing  coherence  between 
water  particle  motion  and  salinity,  sound  velovity,  and  temperature. 
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